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Abstract. - The relation between the autocorrelation C(t,t w ) and the integrated linear re- 
sponse function x(fi^ui) is studied in the context of the large- N model for phase-ordering 
systems subjected to a shear flow. In the high temperature phase T > T c a non-equilibrium 
stationary state is entered which is characterized by a non-trivial fluctuation-dissipation relation 
x(t — t w ) = x(C(t — t w )). For quenches below T c the splitting of the order parameter field into 
two statistically independent components, responsible for the stationary C st (t — t w ) and aging 
C ag (t/t w ) part of the autocorrelation function, can be explicitly exhibited in close analogy with 
the undriven case. In the regime t — t w <C t w the same relation x(t — t w ) — xiP (t — t w )) 
is found between the response and C 3t (t — t iu), as for T > T c . The aging part of x(^>^ ; ) is 
negligible for t m — > oo, as without drive, resulting in a flat x(C) i n the aging regime t — t w 2> t w . 



Understanding non-equilibrium processes is a central topic of modern statistical physics. 
Many systems in nature are observed in a persistent non-equilibrium regime. Among these 
two important categories can be identified according to the asymptotic state towards which 
they evolve. Systems subjected to a change of thermodynamic variables, such as temperature 
or pressure, belong to the first class. In this case, an equilibrium state is approached which, 
however, can be hardly reached if the ergodic time grows in the thermodynamic limit. This 
is usually observed when the initial and final states are separated by a phase transition. 
Typical examples are glassy materials or phase-separating systems quenched below the critical 
temperature. The off-equilibrium dynamics in this case is characterized by the aging property 
of the autocorrelation function Jl| . A different situation, on the other hand, is observed when 
external power is injected into a system with rate 7 by doing work on it, as in the case of stirred 
or sheared materials. In this case a time translational invariant (TTI) stationary regime can 
be entered, depending on the thermodynamic parameters, which is not an equilibrium state. 
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A characterization of non-equilibrium states has been given in different contexts by general- 
izing the fluctuation-dissipation theorem (FDT). In equilibrium the FDT relates the integrated 
linear response function x(t ~ tw) to a perturbation applied at time t w and the autocorrelation 
function C(t — t w ) of the unperturbed system through x(< — tw) = (1/T)[C(0) — C[t — t w )\. 
In stationary states these quantities still depend on the time difference so that \{t — t w ) = 
x[C(t — t w )], where x is a non trivial function that may give informations on the off-equilibrium 
properties. Also in the late stage evolution of aging systems the dependence of the response 
function on the two times generally enters through the autocorrelation, x{t,t w ) — x[C(t,t w )], 
as originally observed by Cugliandolo and Kurchan in the framework of mean field glass 
models. 

In these cases the fluctuation dissipation ratio (FDR) X(C) = —Tdx{C)/dC describes 
deviations from equilibrium, where X(C) = I, and different physical informations can be 
extracted from it. First, T eff = TX' 1 ^) can be interpreted as an effective temperature of 
the system, different from that of the thermal bath, as recognized in turbulence || or glass 
physics ||, [|. Second, the theorem by Franz, Mezard, Parisi and Peliti [[| connects the FDR of 
aging systems to the static properties of the equilibrium state through dX(C) / dC\c=q — P(q), 
where P(q) is the equilibrium probability distribution of the overlaps |6j. Moreover it was 
argued recently [tJ] that aging and driven systems can be related through the off-equilibrium 
generalization of the FDT because X(C) is expected to be the same in the limit of small 
entropy production, namely for t w — * oo in aging systems or 7 — > in driven systems. 

In this Letter we study this subject in the framework of the large-A model for systems 
with non conserved order parameter || under plain shear flow; this is a solvable theory 
which retains the basic features of the equilibrium and off-equilibrium kinetics of several 
non-disordered systems, such as liquid crystals ||, which are relevant for many applications. 
In this context the autocorrelation and the response functions are exactly computed and 
the fluctuation dissipation ratio discussed explicitly. The model is particularly suited for 
a unitary discussion of the FDT in the different non-equilibrium situations discussed above, 
which correspond to dynamical states at different temperatures. In particular, for T larger than 
a critical temperature T c (p/), a stationary state is entered, for long enough times, characterized 
by a non-trivial X(C). By quenching below T c (j), instead, the complex interplay between 
the phase-ordering process and the external drive gives rise to peculiar aging properties. 
Nonetheless, by recognizing that a separation of the degrees of freedom into two statistically 
independent contributions, responsible for the stationary C~(t — t w ) and aging C® 9 (t/t w ) 
part of the autocorrelation, can be explicitly performed, the interpretation of the dynamics is 
clarified by the analogy with the phase-ordering kinetics of the undriven case. 

The model describes the system in terms of an A-component vectorial order parameter field 
<p(r, t) through the Ginzburg- Landau hamiltonian 



where r, u and k are positive constants. When a convective field v is imposed the dynamics 
of the generic component a of a non-conserved order parameter is described by the equation 




(1) 



with 




(2) 



dp, 



'a 



+ V • (<p a v) 




(3) 
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where T is a transport coefficient and the gaussian stochastic variable ff, representing thermal 
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noise, obeys (r] a (r,t)) = and {rj a (r,t)r]p(r' ,t')) = 2T5 a ^TS(r - r')S(t - t'), with < . . . > 
denoting ensemble averages. In the following we will set T = 1. The velocity v is chosen to be 
a planar Couette or shear flow v = jye x , where 7 is the spatially homogeneous shear rate, and 
e x is a unit vector in the flow direction. In the large- N limit summing over infinitely many 
components effectively realizes an ensemble average and the substitution 

1 1 N 

a— 1 

becomes exact JjiJ . Given space-translational invariance the average squared field fluctuation 
S(t) — ((fa) does not depend on the spatial coordinates and (fj) is formally linearized. In 
momentum space it reads 

dJ %^ - ^^Sr = ~ [fc2 + uS[t) ~ rMK t] + ^ t] (5) 

where we have dropped the component index, due to internal symmetry, (rj{k.,t)) = and 
(r](k, t)rj(k',t')) = 2T(27r) D 5(k + k')6(t - t') . The formal solution of Eq.f]) can be obtained 
from the knowledge of the configuration at an arbitrary time to by the method of characteristics 

<P{Kt) =tp[lC(t-t ),to] y-^T e J V~g~{t) e T)[JC(t- z),z]dz 

" t t0 (6) 

where IC(z) = k.+^/zk x e y and g(t) — exp |2 J*[uS(t') — r\dt'\. Eq. (||) allows the computation 
of the structure factor C 7 (k, t) — ((p(k,t)ip(—'k,t)) and the relation 

S(t) = J ^e-^C 7 (k,t), (7) 

where A is a regularization cut-off, provides the self-consistent closure of the model. The 
knowledge of g(i), therefore, amounts to the solution of the theory. 

The existence of critical temperature was shown in [0. For small 7/ A 2 in D = 3 the 
critical temperature behaves as T c {^) — T C (0)[1 -I- ay/ 7A 2 ], where a is a positive constant. In 
general, it can be shown that, in this model, T c (7) is a monotonously increasing function of 
the shear rate. We will consider quenches from an high temperature disordered state to a final 
temperature Tp. In the following we will briefly discuss the most relevant case D — 3. A 
complete analysis, containing also more technical details, will appear elsewhere g(t) can 
be evaluated by Laplace transform techniques, similarly to the case 7 = |l3|, yielding 

r g+ e s °* , for Tp > T c ( 7 ) 
g(t) = { g c , for T F = T C ( 7 ) (8) 

[ 5 _ t~i , forT F <T C ( 7 ), 

where g+,g c ,g- are constants. sq(T) is related to the longitudinal coherence length by sq = 
2£(T 1 ) 2 - These behaviors are reminiscent of the case without shear where the algebraic 
decay of g(t) below the critical temperature is a signature of the power-laws characteristic of 
the phase-ordering process. Then, in the case of a quench from high temperature to Tp below 
T c (7), the effects of shear change the phenomenology of the phase-ordering process but do not 
interrupt it. Moreover, since T c (pf) > T c (0), complete phase-ordering is always observed for 
each Tp < T c (0), regardless of the presence of shear. This is an important difference with 
respect to the case often considered M where the aging process is interrupted by a power 
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input. This is what happens when we shake, for instance, a mixture of water and oil to make 
droplets smaller. Actually the possibility of stabilizing a state with a finite coherence length 
by shearing a phase-ordering system is a debated question |l4|, [l5| which may depend crucially 
on hydrodynamic effects |l6| . 

From the knowledge of g(t) two-time quantities can be computed. 

Tf > T c ( 7 ) 

For quenches to Tp > T c (-f) the autocorrelation function C 7 (t,t w ) = (2ir)~ 2D J d D kd D k' 
(<p(k, t)y>(k', t w )) exp{ — (fc 2 + fc' 2 )/(2A 2 )} is a time translational invariant quantity, C 7 (t,t w ) = 
Gf (i -t w ), for t w » l/s - From (§ one finds 



(8tt) 3 / 2 



l 7 2 z 2 + ^ (t _ tji?) 2 



(*-t m ) 2 
«72 



The next point is the calculation of the linear response function lZ(k, t; k', t') = 
(2tt) d 5((p(k, t)/6h(— k', t') \h(k,t)=Q (with t > t') to a perturbation h(r,t) entering the hamil- 
tonian through — J d D rh(r,t) ■ tp(r,t). From this quantity the integrated response function 
is obtained as x 7 (M«0 = fi w dt' {2tt)- 2D J d D kd D k'Tl{^ t; k', t') exp{-(fc 2 + fc' 2 )/(2A 2 )}. For 
t w 3> 1/sq it can be shown that Xt^i^uO = X 7 *(^ — with 



3^ rt-t,„ /„ i 1 



-3/2 



l {t-t W ) = T -^i .l;r+'- ^' . (10) 



i 7 2 z 2 



Writing the argument of y 7 (£ — ttu) in terms of C?*(f — f w ) through (||) one obtains the 
fluctuation dissipation relation X-y(t ~ tw) = X-yiC^) plotted in Fig. [l| The comparison with 
the case without shear, regulated by equilibrium FDT, shows that for small j(t — t w ) the two 
curves collapse but, for larger time differences the lines with 7^0 bend and lower asymptotic 
values are reached. The behavior of the FDR is shown in the inset of Fig. |l|. Starting from 
X(Q) — 1, this quantity decreases to a minimum and then raises towards X(C) = 1 again 
when C — ► 1. Therefore one finds a different and more complex behavior with respect to the 



two-temperature scenario proposed for driven supercooled fluids |21 
Tf = T c (7) 

For quenches at the critical temperature the asymptotic behavior of the two time quantities 
considered above is the same as for Tp > T c {^), with sq = in (||) and (|l0[). In particular 
Cj(t,t w ) and x-y(t>tw) m the large jt w limit depend only out — t w . 

T F < T c (7) 

The situation is more complex when quenches below T c (^) are considered. In this case, 
using (^), in the large t w limit the autocorrelation function can be separated into two contribu- 
tions, namely C 7 (i, t w ) = C^(t~t w ) + C^ 9 (t/t w ), which depend, respectively, on the difference 
and on the ratio between the two times, as usually found in aging systems. This is the signature 
of a fundamental property of the order parameter field <p(r, t), which decomposes below T c (pf) 
into the sum of two independent stochastic fields <p(r, t) = ip{r, t) + a(r, t) responsible for the 



stationary and aging properties of the kinetics, as discussed in 17 . In systems with a scalar 



order parameter ]18| the physical interpretation of this decomposition goes as follows: given a 



Effects of an external drive 



5 



configuration at time to well inside the late stage scaling regime, one can distinguish degrees 
of freedom i/j(r,i) in the bulk of domains from those cr(r, t) pertaining to the interfaces. The 
first ones are driven by thermal fluctuations and behave locally as in equilibrium; the second 
ones retain memory of the noisy initial condition, and produce the aging behavior. In the 
limit N — > oo the separation of ip(k, t) can be explicitly constructed, identifying in (Q) the two 
contributions 




V(M)=lim / J^e-O^W%[JC(t-z),z]d* (11) 



and 



a(k, t) = lim <p [JC(t t ), to] J ^-e- JT° (12) 

t ^oo W g(t) 

in complete analogy with the case 7 = (l9). For large- N, where topological defects are 
unstable and domains are not well defined, the recognition of the degrees of freedom associated 
to ip(k, t) and cr (k, t) i s not straightforward as for scalar systems. Nevertheless, the possibility 
of the splitting (0,12) of the order parameter indicates that the same fundamental property 
is shared by systems with different N. 

Building the autocorrelation functions of the new fields (3f (t, t w ) = (2it)~ 2D J dPkdPk' 
(ip(k, i)V>(k', t w )) exp{-(fc 2 +£;' 2 )/(2A 2 )} and C*(t, t w ) = (2ir)- 2D J d D kd D k'(a(k, t)a(k.' , t w )) 
exp{— (fc 2 + fc' 2 )/(2A 2 )} it is possible to prove that, for t w /to sufficiently large, they coincide, 
respectively, with C^(t,t w ) and C® 9 (t,t w ). 

Coming back to the two-time observables the TTI part C^{t — t w ) takes the same form as at 
Tp = T c ( 7 ), namely (§) with s = 0. This quantity decays from Cf (0) = M 2 (0) - M 2 {T F ) to 
zero for time differences of order t w . Here M 7 (Tp) = M 7 (0)[1— Tp/T^-y)] 1 ^ 2 is the spontaneous 
magnetization at stationarity. The aging part, instead, is given by 

is _ [2 . , ( _ t\~* 1 

-(1- 



C7( r H2^M 2 (T F )(l + r ) , (l:!) 



i+¥ (i+¥) 

showing that C® 9 stays approximatively constant, C® 9 (t/t w ) ~ M 2 (T F ), for t — t w <C t w and 
decays for large time differences as (t w /t) x , with A = 5/4, as opposed to the undriven case 
where A = 3/4. Computing the response function it is found that the aging part is negligible 
for large t w and xi^^w) behaves as at T c (j), so that ( |Io| ) with sq — is obeyed. 

We turn now to discuss the relation between x 7 (t,t„,) and C 7 (t,t w ). In the case without 
shear Jl|], the equilibrium FDT relates \ st (t ~ t w ) and C st (t - t w ) by T FX st (t - t w ) = 
C st (0) — C st (t — t w ) at every temperature. In the range of times t — t w < L (stationary 
regime), when C st (t - t w ) decays to zero, C ag (t/t w ) ~ M 2 (T F ) and T F x st (t,t w ) grows to its 
saturation value M 2 (0) — M 2 {T F ). Then, plotting T F x(t,t w ) against C(t,t w ) one gets the 
straight line with negative slope on the right of Fig. ||. For larger time differences, in the 
aging regime t — t w t w , x(t^w) is a constant, because the aging part of the response is 
negligible for t w — > 00. In the meanwhile C ag (t/t w ) slowly decreases, and the characteristic 
horizontal line of phase ordering is found. When the drive is switched on, (i — tw) an d 
C^(t — t w ) are not related by equilibrium FDT but, instead, through x 7 (C 7 *) (Fig. [I]). Then 
in the parametric plot of x 7 (t,t„,) vs Cj(t,t w ) (Fig. |^) one finds, in place of the straight line 
with negative slope, the curve x 7 (C 7 *) which corresponds to the simple translation of the one 
plotted in Fig. |l|. For C 7 (t,t w ) < M 2 (T F ) the plot is flat for the same reasons as for 7 = 0. 
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However, due to the bending of x 7 (C*'), the asymptotic value x 7 (oo) is different from the case 
7 = 0. ft is interesting to observe that the relation Tpx(oo) = M 2 (0) — M 2 (Tp) characteristic 
of the undriven case is generalized to Tpx 7 / 2 (oo) = M 7 (0) — M 2 (Tp). Therefore a relation 
analogous to the case 7 = holds between the response function of a system with shear rate 
7/2 and the stationary magnetization in the presence of a shear rate 7 twice as large. This 
parameter-free analytical prediction is amenable of experimental testing. 

In summary, in this Letter we have studied the effects of a shear flow on the fluctuation- 
dissipation relation in the large-iV model for phase-ordering. This approach allows for the 
analytical computation of two times quantities and the explicit identification of the stochastic 
fields -0(r, 4) , a(r,t) responsible for the stationary and aging properties of the kinetics below 
1^(7). The FDT relating the stationary parts of the correlation and response function in the 
undriven case is generalized under shear into a non-trivial relation x 7 '(* — t w ) = X-y(C- st (t — 
t w )), which is obeyed at every temperature. Below T c (^) the lack of time translational 
invariance is testified by the aging contribution C ag (t/t w ) to the correlation function, which 
behaves asymptotically as (t w /t) x , where A = (D + 2)/4 is a generalization of the Fisher-Husc 
exponent |^0| to sheared systems. The aging part of the response, on the other hand, is always 
negligible in the large t w limit, resulting in a flat part of the fluctuation dissipation plot. The 
global picture provided by the model is closely related to the one observed without shear, the 
most notable modification being the different relation between X?Ht — t w ) and C 7 *(t — t w ). 
The comparison with the case 7 = will be discussed further in J12[. 



The applicability of the two-temperature scenario proposed in J21| for supercooled fluids 
to different slowly relaxing systems and, in particular, to the phase-ordering kinetics is a 
debated Jl6| issue which is worth of further studies. According to this scenario an external 
input of energy stabilizes an aging system into a power-dependent stationary state. In the 
small 7 limit the fluctuation dissipation plot of this state is a broken line. This provides the 
identification of two well defined temperatures: The bath temperature Tf and an effective 
temperature of the system. The results of this Letter allows this problem to be addressed in 
the context of the large- iV model. We find that a phase-ordering system is not stabilized by 
shear. This result agrees with Rcf. [|l5| , where a different approximation scheme is developed. A 
similar behavior was found in ]22] ] for the stability of the ferromagnetic phase of the asymmetric 
spin-glass model in the mean- field limit. This feature is at variance with the basic assumption 
of the scenario proposed for supercooled fluids. Furthermore, the behavior of X(C) shown in 
the inset of Fig. [l] shows that a two-temperature interpretation of the fluctuation dissipation 
relation does not hold in this case. 

The relevance of these results to real physical systems where N is finite and hydrodynamic 
effects may play a relevant role is a crucial point that deserves further investigations. 
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C/(t-t w ) 

Fig. 1. - TFX^(t — t w ) is plotted against C^{t — t w ) for a quench to a final temperature 2> = 
(21/20)T C (7), with t w = 10 5 and A = 1. Solid lines correspond to 7 = 1, 10"\ 1CT 2 , bottom up. The 
dashed line, which is almost indistinguishable from the case 7 = 1CP 2 , is the equilibrium behavior 
7 = 0. The inset shows X(C). 
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C y (U w ) 

Fig. 2. - TV Y 7 (*,£,») is plotted against C y (t,t w ) for a quench to a final temperature Tf = T c (j)/2, 
with t w = 10 and A = 1. Solid lines correspond to 7 = 1, 10 _1 , 10 -2 , bottom up. The dashed line is 
the case with 7 = 0. 



